Phenotype-based targeted treatment of SGLT2 inhibitors and GLP-1 receptor agonists in type 2 diabetes.

AIMS/HYPOTHESIS: A precision medicine approach in type 2 diabetes could enhance targeting specific glucose-lowering therapies to individual patients most likely to benefit. We aimed to use the recently developed Bayesian causal forest (BCF) method to develop and validate an individualised treatment selection algorithm for two major type 2 diabetes drug classes, sodium-glucose cotransporter 2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP1-RA). METHODS: We designed a predictive algorithm using BCF to estimate individual-level conditional average treatment effects for 12-month glycaemic outcome (HbA1c) between SGLT2i and GLP1-RA, based on routine clinical features of 46,394 people with type 2 diabetes in primary care in England (Clinical Practice Research Datalink; 27,319 for model development, 19,075 for hold-out validation), with additional external validation in 2252 people with type 2 diabetes from Scotland (SCI-Diabetes [Tayside & Fife]). Differences in glycaemic outcome with GLP1-RA by sex seen in clinical data were replicated in clinical trial data (HARMONY programme: liraglutide [n=389] and albiglutide [n=1682]). As secondary outcomes, we evaluated the impacts of targeting therapy based on glycaemic response on weight change, tolerability and longer-term risk of new-onset microvascular complications, macrovascular complications and adverse kidney events. RESULTS: Model development identified marked heterogeneity in glycaemic response, with 4787 (17.5%) of the development cohort having a predicted HbA1c benefit >3 mmol/mol (>0.3%) with SGLT2i over GLP1-RA and 5551 (20.3%) having a predicted HbA1c benefit >3 mmol/mol with GLP1-RA over SGLT2i. Calibration was good in hold-back validation, and external validation in an independent Scottish dataset identified clear differences in glycaemic outcomes between those predicted to benefit from each therapy. Sex, with women markedly more responsive to GLP1-RA, was identified as a major treatment effect modifier in both the UK observational datasets and in clinical trial data: HARMONY-7 liraglutide (GLP1-RA): 4.4 mmol/mol (95% credible interval [95% CrI] 2.2, 6.3) (0.4% [95% CrI 0.2, 0.6]) greater response in women than men. Targeting the two therapies based on predicted glycaemic response was also associated with improvements in short-term tolerability and long-term risk of new-onset microvascular complications. CONCLUSIONS/INTERPRETATION: Precision medicine approaches can facilitate effective individualised treatment choice between SGLT2i and GLP1-RA therapies, and the use of routinely collected clinical features for treatment selection could support low-cost deployment in many countries.

Dirichlet process mixture models to impute missing predictor data in counterfactual prediction models: an application to predict optimal type 2 diabetes therapy.

BACKGROUND: The handling of missing data is a challenge for inference and regression modelling. A particular challenge is dealing with missing predictor information, particularly when trying to build and make predictions from models for use in clinical practice. METHODS: We utilise a flexible Bayesian approach for handling missing predictor information in regression models. This provides practitioners with full posterior predictive distributions for both the missing predictor information (conditional on the observed predictors) and the outcome-of-interest. We apply this approach to a previously proposed counterfactual treatment selection model for type 2 diabetes second-line therapies. Our approach combines a regression model and a Dirichlet process mixture model (DPMM), where the former defines the treatment selection model, and the latter provides a flexible way to model the joint distribution of the predictors. RESULTS: We show that DPMMs can model complex relationships between predictor variables and can provide powerful means of fitting models to incomplete data (under missing-completely-at-random and missing-at-random assumptions). This framework ensures that the posterior distribution for the parameters and the conditional average treatment effect estimates automatically reflect the additional uncertainties associated with missing data due to the hierarchical model structure. We also demonstrate that in the presence of multiple missing predictors, the DPMM model can be used to explore which variable(s), if collected, could provide the most additional information about the likely outcome. CONCLUSIONS: When developing clinical prediction models, DPMMs offer a flexible way to model complex covariate structures and handle missing predictor information. DPMM-based counterfactual prediction models can also provide additional information to support clinical decision-making, including allowing predictions with appropriate uncertainty to be made for individuals with incomplete predictor data.

Body mass index and inflammation in depression and treatment-resistant depression: a Mendelian randomisation study.

BACKGROUND: Major depressive disorder (MDD) has a significant impact on global burden of disease. Complications in clinical management can occur when response to pharmacological modalities is considered inadequate and symptoms persist (treatment-resistant depression (TRD)). We aim to investigate inflammation, proxied by C-reactive protein (CRP) levels, and body mass index (BMI) as putative causal risk factors for depression and subsequent treatment resistance, leveraging genetic information to avoid confounding via Mendelian randomisation (MR). METHODS: We used the European UK Biobank subcohort ([Formula: see text]), the mental health questionnaire (MHQ) and clinical records. For treatment resistance, a previously curated phenotype based on general practitioner (GP) records and prescription data was employed. We applied univariable and multivariable MR models to genetically predict the exposures and assess their causal contribution to a range of depression outcomes. We used a range of univariable, multivariable and mediation MR models techniques to address our research question with maximum rigour. In addition, we developed a novel statistical procedure to apply pleiotropy-robust multivariable MR to one sample data and employed a Bayesian bootstrap procedure to accurately quantify estimate uncertainty in mediation analysis which outperforms standard approaches in sparse binary outcomes. Given the flexibility of the one-sample design, we evaluated age and sex as moderators of the effects. RESULTS: In univariable MR models, genetically predicted BMI was positively associated with depression outcomes, including MDD ([Formula: see text] ([Formula: see text] CI): 0.133(0.072, 0.205)) and TRD (0.347(0.002, 0.682)), with a larger magnitude in females and with age acting as a moderator of the effect of BMI on severity of depression (0.22(0.050, 0.389)). Multivariable MR analyses suggested an independent causal effect of BMI on TRD not through CRP (0.395(0.004, 0.732)). Our mediation analyses suggested that the effect of CRP on severity of depression was partly mediated by BMI. Individuals with TRD ([Formula: see text]) observationally had higher CRP and BMI compared with individuals with MDD alone and healthy controls. DISCUSSION: Our work supports the assertion that BMI exerts a causal effect on a range of clinical and questionnaire-based depression phenotypes, with the effect being stronger in females and in younger individuals. We show that this effect is independent of inflammation proxied by CRP levels as the effects of CRP do not persist when jointly estimated with BMI. This is consistent with previous evidence suggesting that overweight contributed to depression even in the absence of any metabolic consequences. It appears that BMI exerts an effect on TRD that persists when we account for BMI influencing MDD.

Multi-locus homozygosity promotes actuarial senescence in a wild mammal.

Genome-wide homozygosity, caused for example by inbreeding, is expected to have deleterious effects on survival and/or reproduction. Evolutionary theory predicts that any fitness costs are likely to be detected in late life because natural selection will filter out negative impacts on younger individuals with greater reproductive value. Here we infer associations between multi-locus homozygosity (MLH), sex, disease and age-dependent mortality risks using Bayesian analysis of the life histories of wild European badgers Meles meles in a population naturally infected with Mycobacterium bovis (the causative agent of bovine tuberculosis [bTB]). We find important effects of MLH on all parameters of the Gompertz-Makeham mortality hazard function, but particularly in later life. Our findings confirm the predicted association between genomic homozygosity and actuarial senescence. Increased homozygosity is particularly associated with an earlier onset, and greater rates of actuarial senescence, regardless of sex. The association between homozygosity and actuarial senescence is further amplified among badgers putatively infected with bTB. These results recommend further investigation into the ecological and behavioural processes that result in genome-wide homozygosity, and focused work on whether homozygosity is harmful or beneficial during early life-stages.

Demonstrating multi-country calibration of a tuberculosis model using new history matching and emulation package - hmer.

Infectious disease models are widely used by epidemiologists to improve the understanding of transmission dynamics and disease natural history, and to predict the possible effects of interventions. As the complexity of such models increases, however, it becomes increasingly challenging to robustly calibrate them to empirical data. History matching with emulation is a calibration method that has been successfully applied to such models, but has not been widely used in epidemiology partly due to the lack of available software. To address this issue, we developed a new, user-friendly R package hmer to simply and efficiently perform history matching with emulation. In this paper, we demonstrate the first use of hmer for calibrating a complex deterministic model for the country-level implementation of tuberculosis vaccines to 115 low- and middle-income countries. The model was fit to 9-13 target measures, by varying 19-22 input parameters. Overall, 105 countries were successfully calibrated. Among the remaining countries, hmer visualisation tools, combined with derivative emulation methods, provided strong evidence that the models were misspecified and could not be calibrated to the target ranges. This work shows that hmer can be used to simply and rapidly calibrate a complex model to data from over 100 countries, making it a useful addition to the epidemiologist's calibration tool-kit.

Weak and pleiotropy robust sex-stratified Mendelian randomization in the one sample and two sample settings.

BACKGROUND: Mendelian randomization (MR) leverages genetic data as an instrumental variable to provide estimates for the causal effect of an exposure X on a health outcome Y that is robust to confounding. Unfortunately, horizontal pleiotropy-the direct association of a genetic variant with multiple phenotypes-is highly prevalent and can easily render a genetic variant an invalid instrument. METHODS: Building on existing work, we propose a simple method for leveraging sex-specific genetic associations to perform weak and pleiotropy-robust MR analysis. This is achieved by constructing an MR estimator in which pleiotropy is perfectly removed by cancellation, while placing it within the powerful machinery of the robust adjusted profile score (MR-RAPS) method. Pleiotropy cancellation has the attractive property that it removes heterogeneity and therefore justifies a statistically efficient fixed effects model. We extend the method from the typical two-sample summary-data MR setting to the one-sample setting by adapting the technique of Collider-Correction. Simulation studies and applied examples are used to assess how the sex-stratified MR-RAPS estimator performs against other common approaches. RESULTS: The sex-stratified MR-RAPS method is shown to be robust to pleiotropy even in cases where all genetic variants violated the standard Instrument Strength Independent of Direct Effect assumption. In some cases where the strength of the pleiotropic effect additionally varied by sex (and so perfect cancellation was not achieved), over-dispersed MR-RAPS implementations can still consistently estimate the true causal effect. In applied analyses, we investigate the causal effect of waist-hip ratio (WHR), an important marker of central obesity, on a range of downstream traits. While the conventional approaches suggested paradoxical links between WHR and height and body mass index, the sex-stratified approach obtained a more realistic null effect. Nonzero effects were also detected for systolic and diastolic blood pressure as well as high-density and low-density lipoprotein cholesterol. DISCUSSION: We provide a simple but attractive method for weak and pleiotropy robust causal estimation of sexually dimorphic traits on downstream outcomes, by combining several existing approaches in a novel fashion.

Islet Autoantibody Level Distribution in Type 1 Diabetes and Their Association With Genetic and Clinical Characteristics.

CONTEXT: The importance of the autoantibody level at diagnosis of type 1 diabetes (T1D) is not clear. OBJECTIVE: We aimed to assess the association of glutamate decarboxylase (GADA), islet antigen-2 (IA-2A), and zinc transporter 8 (ZnT8A) autoantibody levels with clinical and genetic characteristics at diagnosis of T1D. METHODS: We conducted a prospective, cross-sectional study. GADA, IA-2A, and ZnT8A were measured in 1644 individuals with T1D at diagnosis using radiobinding assays. Associations between autoantibody levels and the clinical and genetic characteristics for individuals were assessed in those positive for these autoantibodies. We performed replication in an independent cohort of 449 people with T1D. RESULTS: GADA and IA-2A levels exhibited a bimodal distribution at diagnosis. High GADA level was associated with older age at diagnosis (median 27 years vs 19 years, P = 9 × 10-17), female sex (52% vs 37%, P = 1 × 10-8), other autoimmune diseases (13% vs 6%, P = 3 × 10-6), and HLA-DR3-DQ2 (58% vs 51%, P = .006). High IA-2A level was associated with younger age of diagnosis (median 17 years vs 23 years, P = 3 × 10-7), HLA-DR4-DQ8 (66% vs 50%, P = 1 × 10-6), and ZnT8A positivity (77% vs 52%, P = 1 × 10-15). We replicated our findings in an independent cohort of 449 people with T1D where autoantibodies were measured using enzyme-linked immunosorbent assays. CONCLUSION: Islet autoantibody levels provide additional information over positivity in T1D at diagnosis. Bimodality of GADA and IA-2A autoantibody levels highlights the novel aspect of heterogeneity of T1D. This may have implications for T1D prediction, treatment, and pathogenesis.

Estimating cetacean population trends from static acoustic monitoring data using Paired Year Ratio Assessment (PYRA).

The cetacean conservationist is often faced with evaluating population trends from abundance data that are either sparse or recorded at different times in different years. The presence of diel or seasonal patterns in the data together with unplanned gaps is often problematic. Such data are typical of those obtained from static acoustic monitoring. We present a simple and transparent non-parametric trend evaluation method, 'Paired Year Ratio Assessment (PYRA)' that uses only whole days of data wherever they are present in each of successive pairs of periods of 365 days. We provide a quantitative comparison of the performance of PYRA with traditional generalised additive models (GAMS) and nonparametric randomisation tests that require a greater level of skill and experience for both application and interpretation. We conclude that PYRA is a powerful tool, particularly in the context of identifying population trends which is often the main aim of conservation-targeted acoustic monitoring.

Challenges in estimation, uncertainty quantification and elicitation for pandemic modelling.

The estimation of parameters and model structure for informing infectious disease response has become a focal point of the recent pandemic. However, it has also highlighted a plethora of challenges remaining in the fast and robust extraction of information using data and models to help inform policy. In this paper, we identify and discuss four broad challenges in the estimation paradigm relating to infectious disease modelling, namely the Uncertainty Quantification framework, data challenges in estimation, model-based inference and prediction, and expert judgement. We also postulate priorities in estimation methodology to facilitate preparation for future pandemics.

Spatial and temporal variation in proximity networks of commercial dairy cattle in Great Britain.

The nature of contacts between hosts can be important in facilitating or impeding the spread of pathogens within a population. Networks constructed from contacts between hosts allow examination of how individual variation might influence the spread of infections. Studying the contact networks of livestock species managed under different conditions can additionally provide insight into their influence on these contact structures. We collected high-resolution proximity and GPS location data from nine groups of domestic cattle (mean group size = 85) in seven dairy herds employing a range of grazing and housing regimes. Networks were constructed from cattle contacts (defined by proximity) aggregated by different temporal windows (2 h, 24 h, and approximately 1 week) and by location within the farm. Networks of contacts aggregated over the whole study were highly saturated but dividing contacts by space and time revealed substantial variation in cattle interactions. Cows showed statistically significant variation in the frequency of their contacts and in the number of cows with which they were in contact. When cows were in buildings, compared to being on pasture, contact durations were longer and cows contacted more other cows. A small number of cows showed evidence of consistent relationships but the majority of cattle did not. In one group where management allowed free access to all farm areas, cows showed asynchronous space use and, while at pasture, contacted fewer other cows and showed substantially greater between-individual variation in contacts than other groups. We highlight the degree to which variations in management (e.g. grazing access, milking routine) substantially alter cattle contact patterns, with potentially major implications for infection transmission and social interactions. In particular, where individual cows have free choice of their environment, the resulting contact networks may have a less-risky structure that could reduce the likelihood of direct transmission of infections.

Long-term effects of antibiotic treatments on honeybee colony fitness: A modelling approach.

Gut microbiome disequilibrium is increasingly implicated in host fitness reductions, including for the economically important and disease-challenged western honey bee Apis mellifera. In laboratory experiments, the antibiotic tetracycline, which is used to prevent American Foulbrood Disease in countries including the US, elevates honey bee mortality by disturbing the microbiome. It is unclear, however, how elevated individual mortality affects colony-level fitness.We used an agent-based model (BEEHAVE) and empirical data to assess colony-level effects of antibiotic-induced worker bee mortality, by measuring colony size. We investigated the relationship between the duration that the antibiotic-induced mortality probability is imposed for and colony size.We found that when simulating antibiotic-induced mortality of worker bees from just 60 days per year, up to a permanent effect, the colony is reduced such that tetracycline treatment would not meet the European Food Safety Authority's (EFSA) honey bee protection goals. When antibiotic mortality was imposed for the hypothetical minimal exposure time, which assumes that antibiotics only impact the bee's fitness during the recommended treatment period of 15 days in both spring and autumn, the colony fitness reduction was only marginally under the EFSA's threshold. Synthesis and Applications. Modelling colony-level impacts of antibiotic treatment shows that individual honey bee worker mortality can lead to colony mortality. To assess the full impact, the persistence of antibiotic-induced mortality in honey bees must be determined experimentally, in vivo. We caution that as the domestication of new insect species increases, maintaining healthy gut microbiomes is of paramount importance to insect health and commercial productivity. The recommendation from this work is to limit prophylactic use of antibiotics and to not exceed recommended treatment strategies for domesticated insects. This is especially important for highly social insects as excess antibiotic use will likely decrease colony growth and an increase in colony mortality.

Characterization of potential superspreader farms for bovine tuberculosis: A review.

BACKGROUND: Variation in host attributes that influence their contact rates and infectiousness can lead some individuals to make disproportionate contributions to the spread of infections. Understanding the roles of such 'superspreaders' can be crucial in deciding where to direct disease surveillance and controls to greatest effect. In the epidemiology of bovine tuberculosis (bTB) in Great Britain, it has been suggested that a minority of cattle farms or herds might make disproportionate contributions to the spread of Mycobacterium bovis, and hence might be considered 'superspreader farms'. OBJECTIVES AND METHODS: We review the literature to identify the characteristics of farms that have the potential to contribute to exceptional values in the three main components of the farm reproductive number - Rf : contact rate, infectiousness and duration of infectiousness, and therefore might characterize potential superspreader farms for bovine tuberculosis in Great Britain. RESULTS: Farms exhibit marked heterogeneity in contact rates arising from between-farm trading of cattle. A minority of farms act as trading hubs that greatly augment connections within cattle trading networks. Herd infectiousness might be increased by high within-herd transmission or the presence of supershedding individuals, or infectiousness might be prolonged due to undetected infections or by repeated local transmission, via wildlife or fomites. CONCLUSIONS: Targeting control methods on putative superspreader farms might yield disproportionate benefits in controlling endemic bovine tuberculosis in Great Britain. However, real-time identification of any such farms, and integration of controls with industry practices, present analytical, operational and policy challenges.

Key questions for modelling COVID-19 exit strategies.

Combinations of intense non-pharmaceutical interventions (lockdowns) were introduced worldwide to reduce SARS-CoV-2 transmission. Many governments have begun to implement exit strategies that relax restrictions while attempting to control the risk of a surge in cases. Mathematical modelling has played a central role in guiding interventions, but the challenge of designing optimal exit strategies in the face of ongoing transmission is unprecedented. Here, we report discussions from the Isaac Newton Institute 'Models for an exit strategy' workshop (11-15 May 2020). A diverse community of modellers who are providing evidence to governments worldwide were asked to identify the main questions that, if answered, would allow for more accurate predictions of the effects of different exit strategies. Based on these questions, we propose a roadmap to facilitate the development of reliable models to guide exit strategies. This roadmap requires a global collaborative effort from the scientific community and policymakers, and has three parts: (i) improve estimation of key epidemiological parameters; (ii) understand sources of heterogeneity in populations; and (iii) focus on requirements for data collection, particularly in low-to-middle-income countries. This will provide important information for planning exit strategies that balance socio-economic benefits with public health.

Effects of trading networks on the risk of bovine tuberculosis incidents on cattle farms in Great Britain.

Trading animals between farms and via markets can provide a conduit for spread of infections. By studying trading networks, we might better understand the dynamics of livestock diseases. We constructed ingoing contact chains of cattle farms in Great Britain that were linked by trading, to elucidate potential pathways for the transmission of infection and to evaluate their effect on the risk of a farm experiencing a bovine tuberculosis (bTB) incident. Our findings are consistent with variation in bTB risk associated with region, herd size, disease risk area and history of previous bTB incidents on the root farm and nearby farms. However, we also identified effects of both direct and indirect trading patterns, such that connections to more farms in the England High-Risk Area up to three movements away from the root farm increased the odds of a bTB incident, while connections with more farms in the England Low-Risk Area up to eight movements away decreased the odds. Relative to other risk factors for bTB, trading behaviours are arguably more amenable to change, and consideration of risks associated with indirect trading, as well direct trading, might therefore offer an additional approach to bTB control in Great Britain.

Apparent absence of Batrachochytrium salamandrivorans in wild urodeles in the United Kingdom.

Whether an infectious disease threat to wildlife arises from pathogen introduction or the increased incidence of an already-present agent informs mitigation policy and actions. The prior absence of a pathogen can be difficult to establish, particularly in free-living wildlife. Subsequent to the epidemic emergence of the fungus, Batrachochytrium salamandrivorans (Bsal), in mainland Europe in 2010 and prior to its detection in captive amphibians in the United Kingdom (UK), we tested archived skin swabs using a Bsal-specific qPCR. These samples had been collected in 2011 from 2409 wild newts from ponds across the UK. All swabs were negative for Bsal. Bayesian hierarchical modelling suggests that Bsal was absent from, or present at very low levels in, these ponds at the time of sampling. Additionally, surveillance of newt mortality incidents, 2013-2017, failed to detect Bsal. As this pathogen has been shown to be widespread in British captive amphibian collections, there is an urgent need to raise awareness of the importance of effective biosecurity measures, especially amongst people with captive amphibians, to help minimise the risk of Bsal spreading to the wild. Continued and heightened wild amphibian disease surveillance is a priority to provide an early warning system for potential incursion events.

Contact chains of cattle farms in Great Britain.

Network analyses can assist in predicting the course of epidemics. Time-directed paths or 'contact chains' provide a measure of host-connectedness across specified timeframes, and so represent potential pathways for spread of infections with different epidemiological characteristics. We analysed networks and contact chains of cattle farms in Great Britain using Cattle Tracing System data from 2001 to 2015. We focused on the potential for between-farm transmission of bovine tuberculosis, a chronic infection with potential for hidden spread through the network. Networks were characterized by scale-free type properties, where individual farms were found to be influential 'hubs' in the network. We found a markedly bimodal distribution of farms with either small or very large ingoing and outgoing contact chains (ICCs and OCCs). As a result of their cattle purchases within 12-month periods, 47% of British farms were connected by ICCs to more than 1000 other farms and 16% were connected to more than 10 000 other farms. As a result of their cattle sales within 12-month periods, 66% of farms had OCCs that reached more than 1000 other farms and 15% reached more than 10 000 other farms. Over 19 000 farms had both ICCs and OCCs reaching more than 10 000 farms for two or more years. While farms with more contacts in their ICCs or OCCs might play an important role in disease spread, farms with extensive ICCs and OCCs might be particularly important by being at higher risk of both acquiring and disseminating infections.

Risk factors and variations in detection of new bovine tuberculosis breakdowns via slaughterhouse surveillance in Great Britain.

Slaughterhouse surveillance through post-mortem meat inspection provides an important mechanism for detecting bovine tuberculosis (bTB) infections in cattle herds in Great Britain (GB), complementary to the live animal skin test based programme. We explore patterns in the numbers of herd breakdowns detected through slaughterhouse surveillance and develop a Bayesian hierarchical regression model to assess the associations of animal-level factors with the odds of an infected animal being detected in the slaughterhouse, allowing us to highlight slaughterhouses that show atypical patterns of detection. The analyses demonstrate that the numbers and proportions of breakdowns detected in slaughterhouses increased in GB over the period of study (1998-2013). The odds of an animal being a slaughterhouse case was strongly associated with the region of the country that the animal spent most of its life, with animals living in high-frequency testing areas of England having on average 21 times the odds of detection compared to animals living in Scotland. There was also a strong effect of age, with animals slaughtered at > 60 months of age having 5.3 times the odds of detection compared to animals slaughtered between 0-18 months of age. Smaller effects were observed for cattle having spent time on farms with a history of bTB, quarter of the year that the animal was slaughtered, movement and test history. Over-and-above these risks, the odds of detection increased by a factor of 1.1 for each year of the study. After adjustment for these variables, there were additional variations in risk between slaughterhouses and breed. Our framework has been adopted into the routine annual surveillance reporting carried out by the Animal Plant Health Agency and may be used to target more detailed investigation of meat inspection practices.

Choice of time horizon critical in estimating costs and effects of changes to HIV programmes.

BACKGROUND: Uganda changed its antiretroviral therapy guidelines in 2014, increasing the CD4 threshold for antiretroviral therapy initiation from 350 cells/µl to 500 cells/µl. We investigate what effect this change in policy is likely to have on HIV incidence, morbidity, and programme costs, and estimate the cost-effectiveness of the change over different time horizons. METHODS: We used a complex individual-based model of HIV transmission and antiretroviral therapy scale-up in Uganda. 100 model fits were generated by fitting the model to 51 demographic, sexual behaviour, and epidemiological calibration targets, varying 96 input parameters, using history matching with model emulation. An additional 19 cost and disability weight parameters were varied during the analysis of the model results. For each model fit, the model was run to 2030, with and without the change in threshold to 500 cells/µl. RESULTS: The change in threshold led to a 9.7% (90% plausible range: 4.3%-15.0%) reduction in incidence in 2030, and averted 278,944 (118,452-502,790) DALYs, at a total cost of $28M (-$142M to +$195M). The cost per disability adjusted life year (DALY) averted fell over time, from $3238 (-$125 to +$29,969) in 2014 to $100 (-$499 to +$785) in 2030. The change in threshold was cost-effective (cost <3×Uganda's per capita GDP per DALY averted) by 2018, and highly cost-effective (cost

Improving ART programme retention and viral suppression are key to maximising impact of treatment as prevention - a modelling study.

BACKGROUND: UNAIDS calls for fewer than 500,000 new HIV infections/year by 2020, with treatment-as-prevention being a key part of their strategy for achieving the target. A better understanding of the contribution to transmission of people at different stages of the care pathway can help focus intervention services at populations where they may have the greatest effect. We investigate this using Uganda as a case study. METHODS: An individual-based HIV/ART model was fitted using history matching. 100 model fits were generated to account for uncertainties in sexual behaviour, HIV epidemiology, and ART coverage up to 2015 in Uganda. A number of different ART scale-up intervention scenarios were simulated between 2016 and 2030. The incidence and proportion of transmission over time from people with primary infection, post-primary ART-naïve infection, and people currently or previously on ART was calculated. RESULTS: In all scenarios, the proportion of transmission by ART-naïve people decreases, from 70% (61%-79%) in 2015 to between 23% (15%-40%) and 47% (35%-61%) in 2030. The proportion of transmission by people on ART increases from 7.8% (3.5%-13%) to between 14% (7.0%-24%) and 38% (21%-55%). The proportion of transmission by ART dropouts increases from 22% (15%-33%) to between 31% (23%-43%) and 56% (43%-70%). CONCLUSIONS: People who are currently or previously on ART are likely to play an increasingly large role in transmission as ART coverage increases in Uganda. Improving retention on ART, and ensuring that people on ART remain virally suppressed, will be key in reducing HIV incidence in Uganda.

Universal test, treat, and keep: improving ART retention is key in cost-effective HIV control in Uganda.

BACKGROUND: With ambitious new UNAIDS targets to end AIDS by 2030, and new WHO treatment guidelines, there is increased interest in the best way to scale-up ART coverage. We investigate the cost-effectiveness of various ART scale-up options in Uganda. METHODS: Individual-based HIV/ART model of Uganda, calibrated using history matching. 22 ART scale-up strategies were simulated from 2016 to 2030, comprising different combinations of six single interventions (1. increased HIV testing rates, 2. no CD4 threshold for ART initiation, 3. improved ART retention, 4. increased ART restart rates, 5. improved linkage to care, 6. improved pre-ART care). The incremental net monetary benefit (NMB) of each intervention was calculated, for a wide range of different willingness/ability to pay (WTP) per DALY averted (health-service perspective, 3% discount rate). RESULTS: For all WTP thresholds above $210, interventions including removing the CD4 threshold were likely to be most cost-effective. At a WTP of $715 (1 × per-capita-GDP) interventions to improve linkage to and retention/re-enrolment in HIV care were highly likely to be more cost-effective than interventions to increase rates of HIV testing. At higher WTP (> ~ $1690), the most cost-effective option was 'Universal Test, Treat, and Keep' (UTTK), which combines interventions 1-5 detailed above. CONCLUSIONS: Our results support new WHO guidelines to remove the CD4 threshold for ART initiation in Uganda. With additional resources, this could be supplemented with interventions aimed at improving linkage to and/or retention in HIV care. To achieve the greatest reductions in HIV incidence, a UTTK policy should be implemented.